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ABSTRACT 

We present a statistical analysis of the properties of galaxies containing ultraluminous X-ray objects 
(ULXs). Our primary goal is to establish the fraction of galaxies containing a ULX as a function of 
ULX luminosity. Our sample is based on ROSAT HRI observations of galaxies. We find that ~ 12% 
of galaxies contain at least one ULX with Lx > 10 39 erg s _1 and ~ 1% of galaxies contain at least one 
ULX with L x > 10 40 erg s" 1 . These ULX frequencies are lower limits since ROSAT HRI observation 
would miss absorbed ULXs (i.e., with Nu > 10 21 cm^ 2 ) and those within ~ 10" of the nucleus (due 
to the positional error circle of the ROSAT HRI). The Hubble type distribution of galaxies with a 
ULX differs significantly from the distribution of types for nearby RC3 galaxies, but does not differ 
significantly from the galaxy type distribution of galaxies observed by the HRI in general. We find 
no increase in the mean FIR luminosity or FIR / K band luminosity ratio for galaxies with a ULX 
relative to galaxies observed by the HRI in general, however this result is also most likely biased by 
the soft bandpass of the HRI and the relatively low number of high SFR galaxies observed by the HRI 
with enough sensitivity to detect a ULX. 
Subject headings: 



1. INTRODUCTION 

ROSAT HRI observations showed that ultraluminous 
X-ray objects (ULXs, a.k.a. IXOs), which were pre - 
viously detected with Einstein (see iFabbianol [l989), 
are compact X-ray sources, are quite common in 
the local Universe, and are typically offset from the 
optical galactic nucleus ijColbert fc Mushotzkvl 119991: 
iRoberts fc War wield 1200(1 . Here we define ULXs as X- 
ray sources in galaxies that are point-like, extra- nuclear, 
and have X-ray luminosities in excess of 10 39 ergs s _1 . 
The re is considerable debate as to nature of these sources 
(see Ivan der Marelll2003l iMiller fc Colbert! 1200.1 for re- 
views), with popular explanations including anisotropic 
emission from stellar-mass X-ray binaries and intermedi- 
ate mass black holes. A constraint on any explanation for 
these sources would be the frequency at which they occur, 
i.e., what fraction of galaxies harbor a ULX as a function 
of ULX luminosity. For example, this frequency would 
have implications for the ty pical relativistic boosting fac- 
tor fo r jet models of ULXs ijKoerding. Falcke. fc Markofj 
2002). In this paper we examine the frequency of ULX 
occurrence as a function of ULX luminosity based on th e 
HRI ULX catalog published in iColbert fc Ptald (|2002|l 
While HRI observations suffer from significant biases 
against the detection of ULXs (as discussed below), we 
note that the relatively large field of view of the HRI 
resulted in a large number of galaxies having been ob- 
served, and of course large area surveys are crucial for 
establishing the bright end of a luminosity distribution. 

2. METHODOLOGY 

In addition to IColbert fc Ptald l|2002ft . we also 
mad e use of a subset of th e RC3 galaxy cata- 
log (|de Vaucouleurs et alJ 11991(1 where the distance 
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is determined, most often based on iTullvl l)1988|) . 
IColbert fc Ptald l(2002fl chose a recessional velocity cut- 
off of 5000 km s _1 . In this paper we assume the cosmol- 
ogy H — 75 km s -1 Mpc" 1 and q = 0.5 (for consistency 
with Tully 1988), and therefore this velocity corresponds 
to 66.7 Mpc. This resulted in a galaxy sample contain- 
ing 9452 galaxies. We next determined which of these 
galaxies had "coverage" by computing the overlap of the 
RC3 D25 ellipse and the HRI FOV for each observation 
that was used in IColbert fc Ptakl l)2002jh We considered 
a galaxy as having coverage if > 90% of the galaxy was 
in the FOV. We found that 766 galaxies had coverage, 
or - 8% of the RC3 subset. 

In order to properly assess the frequency of ULXs, we 
need to assess how many galaxies were observed by the 
HRI with sufficient sensitivity to detect a ULX. We com- 
puted the flux limit at the position of each RC3 galaxy 
with coverage by computing the background level in the 
point sourc e region size. The poin t source region used 
by XAssist ijPtak fc Griffithsll2003D is a circle with a ra- 
dius given by 4[exp(0/9) + 1] arcseconds, where 9 is the 
off-axis angle of the source in arcminutes. The total num- 
ber of background counts, b, in this area was then de- 
termined using the median background counts per pixel 
(also computed as a part of the XAssist processing of 
each field). The 3a detection threshold count rate is 
given by 3y/b/t, where t is the exposure time, and the 
limiting flux is then 3Vbf conv v(8)/t where v(ff) gives the 
vignetting correction and f CO nv is the on-axis count rate 
to flux conversion factor. f conv was computed using the 
FTOOL PIMMS, assuming a power-law slope of 1 . 7 and 
the Galactic column density l|Dickev fc Lockmanlll99fl 
using the FTOOL "nh"). We computed the limiting flux 
in both the 0.5-2.0 keV bandpass F nm,0. 5-2.0 kcV and the 
2-10 keV bandpass F i im ,2-io kov- We used F um.a-io kcV 
to compute the limit ing luminosity s i nce th e 2-10 keV 
bandpass was used in IColbert feTtakl ((2002), for consis- 
tency with current Chandra and XMM-Newton surveys. 
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Any ULX could then have been detected in a galaxy 
when L Um,2-io kcV is < 10 39 erg s" 1 , while more gener- 
ally L u ro ,2-io kev is > 10 39 erg s -1 . 

Finally, it is necessary to determine how many ULXs 
are likely to be background sources . For th i s pur - 
pose we used the logN-logS given in iHasingerl (l998), 
which we note is consistent with a simple power-law 
logiV(> S) ~ lOCKSLur 1 - 4 sources deg" 2 , where 5_ i4 
is the 0.5-2.0 keV flux in units of 10 -14 ergs cm" 2 s" 1 . 
The number of expected background sources is then 
n xrb = 100( F " m ' 1 ° 5 _- 1 2 4 keV )- 1A irr 2 , wh ere r is the ma- 
j or ax is radius of each galaxy in degrees. IColbert fc Ptakl 
( 2002) used r = D25 (the RC3 major axis diameter) how- 
ever here we also consider r = 0.5Z?25- Also note that 
when L n m ,2-io kcV is < 10 39 erg s -1 , the relevant num- 
ber of background sources is the number of source whose 
apparent luminosity exceeds 10 39 erg s _1 , and since we are 
assuming a simple power-law logN-logS this simply en- 
tails decreasing the total number of background sources 
expected by the factor (L n m ,2-io kcv/10 39 erg s" 1 ) 1 - 4 . 

3. RESULTS 

3.1. ULX Frequency 

In Table^we show the number of galaxies {NL Um <Lx) 
with HRI coverage and a limiting 2-10 keV luminosity 
less than or equal to Lx ranging from 10 39 erg s to 
10 41 erg s _1 (i.e., the number of galaxies where a ULX 
as bright as Lx or brighter could have been detected). 
We also show the number of galaxies with at least one 
ULX (N w uLx)i and the number of galaxies with at least 
one ULX whose luminosity exceeds the limiting lumi- 
nosity (N LuLX > Lx ). This latter quantity, N LuLX > Lx , 
is of interest since it represents the number of galaxies 
with a ULX at least as bright as the list Lx, and we 
then computed the resultant number of galaxies after 
subtracting the expected number of background sources 
(N ne t,L UL x>Lx > i- e -i the number of galaxies where the 
1-CT lower limit on the number of ULXs is greater than 
0.; the lower limit was calculated using the mehology de- 
scribed in Kraft, Burrows, & Nousek (1991), with the 
background being the expected number of background 
sources). Finally we list the values N w ulx, ^l ulx >l X i 
and N netyLuLX > Lx using both r = D 25 and r = 0.5L> 2 5- 
No te that in the latter cas e the total number of ULXs in 
the IColbert &: Ptakl l)2002|) sample is reduced from 87 to 
49. 

In Table^shows that a large fraction of ULXs detected 
at radii larger 0.5£>25 are consistent with background 
sources, i.e., the net number of galaxies with ULXs is of- 
ten larger when O.5D25 is chosen as the limiting radius. 
Also, in some cases, the number of ULXs detected within 
the 0.5Z?25 radius is above the background level when r = 
0.5Z?25 but is not above the background level when r = 
D25 (i-e., since in the latter case the expected number of 
background sources is a factor of four larger). In Figure 
Q] we plot the fraction of galaxies harboring a ULX as a 
function of ULX luminosity {N nettLuLX > Lx /N Lli7it < Lx ) 
for both r = D25 and r = 0.5Z?25- Here the errors were 
computed according to the me thodology described in 
iKraft. Burrows, fc Nouse"kl l)1991f ). with the background 
taken as the number of galaxies at each luminosity that 
are consistent with all ULXs being background sources 
(i.e., N L 

ULX>Lx ulx>Lx)- The choice of effec- 
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Fig. 1. — The frequency of ULX occurrence in galaxies as a 
function of ULX luminosity. The frequency of ULXs is plotted 
using effective galactic radii of D25 (blue diamonds with dashed 
error bars) and O.5D25 (red crosses with solid error bars). We also 
plot the frequency observed after limiting the sample to only spirals 
(with Hubble type >= and radius = 0.5-D25; black filled circles 
with dotted error bars). 



tive radius for the galaxies does not significantly im- 
pact the estimate of ULX frequency once background 
sources are taken into account. The frequency of ULXs 
observed by the HRI evidently is ~ 12% and ~ 1% of 
galaxies harboring a ULX with Lx >= 10 39 erg s _1 and 
Lx >= 10 40 erg s _1 , respectively. We also plot the 
galaxy fractions for spirals only (i.e., Hubble type in RC3 
> 0). The Hubble type distribution of galaxies with and 
without a ULX is discussed below, however we note that 
the fraction of galaxies harboring a ULX is slightly higher 
in spirals. 

3.2. Hubble Type Distribution 

In Figure |3 we show the Hubble type distribution of 
the RC3 subsample (discussed in §2 ), the Hubble type 
distribution of galaxies that had HRI coverage (hereafter 
HRI galaxies) and the HRI galaxies that contain at least 
one ULX. The ULX galaxy type distribution is shown 
both before and after background subtraction, and here 
we only consider r < 0.5Z?25 ULXs. There is a slight ex- 
cess of elliptical ULX galaxies relative to the RC3 sam- 
ple however the larger excess in the HRI galaxy sam- 
ple suggests that this bias due in part to a bias toward 
early type galaxies in the HRI observing patterns rela- 
tive to RC3 galaxies in general. The most conspicuous 
feature of the ULX galaxy distribution is the excess of 
T=5-7 galaxies, corresponding to Sc and Scd galaxies. 
However, in general the ULX galaxy type distributions 
are not significantly different statistically from the HRI 
RC3 galaxy type distributions, with the x 2 tests giving 
X 2 /dof ~ 15/8, which can only be rejected at the 95% 
level. 

We also investigated the possible connection between 
star-formation rate (SFR) and ULX frequency by tabu- 
lating the far-infrared luminosities (a well-known SFR 
indicator; Kennicutt 1998) of the galaxy samples dis- 
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Table 1. Galaxy Number Counts 



Lx 
(1) 


N L llpT <L x 
(2) 


N W ULX 

(3) 


N L 


r = D25 
(4) 


Nnet,Lu , X >L X 

(5) 


N W ULX 

(6) 


r 

N L UL 
C 


= 0.5£>25 

X>L X ^net,LuLx> L X 

') (8) 
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228 


45 


45 




25 (8.7%) 


34 


34 


28 (12.3%) 
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316 


49 


38 




23 (7.3%) 
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23 (7.3%) 
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13 (2.9%) 


38 
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14 (3.1%) 
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54 
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38 
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6 (1.1%) 
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1 (0.1%) 


38 





(0.0%) 
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762 
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38 





(0.0%) 



Note. — (1) minimum ULX 2-10 keV luminosity in 10 erg s . (2) JVx, <L X — number of galaxies with a limiting luminosity Ln m less 
than or equal to L x . (3) and (6) N w ulx — number of galaxies with Ln m < L x and a ULX, for r — D25 and ro. 5.D25 ' rcs P cc ti vc ly- (4) and (7) 
N Ljj Lx > L x — number of galaxies with a Lu m < L x and a ULX with luminosity Lulx > Lx- (5) and (8) N net L ULX >L x — number of galaxies 
with a Lnjn < Lx and a net number of ULXs with luminosity Lulx > Lx > after subtracting the expected number of background sources, 
with the percentage of galaxies (i.e., N net l tjlx >l x^ L l' <L x^ S^ ven m parentheses. 
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Fig. 2. — (left) The Hubble type distribution of nearby galaxies in the RC3 catalog (black solid line), the galaxies with HRI coverage 
(red dashed line), the galaxies with a ULX (blue dot-dashed line) and the galaxies with a ULX after excluding galaxies with at least one 
"net" ULX (i.e., after subtracting the expected number of background sources; blue dotted line), (right) The Hubble type histograms 
for galaxies with HRI cove r age (r ed, solid error bars) and with a "net" ULX (blue, dotted error bars) shown with errors computed using 
Kraft. Burrows. & Nousck 1 1991). The galaxy types T correspond to: ellipticals have T =-5 to T=-4, lenticulars have T=-3 to T=-l, 
spirals have T=0 (SO) to T=9, with odd numbers from 1 to 9 representing Sa,b,c,d and even numbers representing Sab, be, cd, dm. Irregulars 
have T=10-ll and we assigned 16 to galaxies with T>11. Both plots were created using only ULXs within r = 0.5-D25- 



cussed above. These distributions are shown in Fig- 
ure where we used NED IRAS fluxes densities at 60 
and 100 ,um to compute F FIR = 1.26 x 10~ 23 [2.58 x 
10 _12 /60/-m + 10~ 12 /iooA I m]- The peak of the distribu- 
tions with a ULX is slightly skewed toward higher FIR 
(i.e., near log L fir = 43.5) however the means of the dis- 
tributions do not differ significantly (all are in the range 
log Lfir = 42.7-43.0). The HRI galaxy FIR luminos- 
ity distribution differs from the "net" ULX galaxy FIR 
luminosity distribution at less than the 90% level based 
on the x 2 test. Restricting the sample to spiral galax- 
ies (T >= 0) shifted the means by only logLpm ~ 0.1, 
reflecting the predominance of spirals in these samples. 

Similarly we plot the distribution of Lfir/Lk where 
Lk is the K band luminosity computed using the 
2mass Extended S ource catalog K band magnitude (see 
IJarrett et alJl2003|) and assuming a zero-point flux den- 



sity of 6.7 x 10~ 21 erg cm" 2 s^Hz' 1 l|Zombeckl ll990) 3 . 
The K band luminosity is a proxy for stellar mass since 
the variation in mass-to-light ratio for K band is only a 
facto r of ~ 2 for various galaxy types (Be ll fc de Jond 
1200 1|) . These distributions are shown in Figure 0] and 
again we do not see any significant differences between 
galaxies with ULXs and galaxies observed by the HRI in 
general. 

4. SUMMARY & CONCLUSIONS 

We have assessed the frequency of ULX occurrence 
in galaxies as a function of ULX luminosity, including 
proper assessment of the contamination from background 
sources. We have found that assuming an effective radius 
of D25 tends to result in a significant contamination from 

3 see also 

://www. ipac.caltech.edu/2mass/releases/allsky/doc/sec6_4a.html 
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Fig. 3. — The FIR luminosity distribution of nearby galaxies with IRAS detections at 60 fim and 100 fim (black solid line), the galaxies 
with HRI coverage (red dashed line), and the galaxies with a ULX before (blue dot-dashed line) and after (blue dotted line) background 
subtraction. The left panel shows the results without selecting on galaxy type and the right panel shows the results when selecting only 
spiral galaxies (Hubble type > 0). 
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Fig. 4. — The FIR/K band luminosity distributions for nearby galaxies with IRAS detections at 60 fim and 100 fim and 2mass extended 
source catalog K band detections (black solid line). Other curves are as in Figure |3 The left panel shows the results without selecting on 
galaxy type and the right panel shows the results when selecting only spiral galaxies (Hubble type > 0). 



background sources. We found frequencies of ~ 12% and 
~ 1% for ULX luminosities of at least 10 39 erg s _1 and 
10 40 erg s _1 . The Hubble type distribution of galaxies 
containing at least one ULX is significantly different from 
those that do not contain a ULX in general, however it 
is not significantly different from galaxies observed by 
the HRI. This implies that any bias in the Hubble type 
distribution of galaxies observed by HRI with ULXs is 
mostly due to the o bserv ing patterns of the HRI. Re- 
cently llrwin et al. I l|2004|) reported on an analysis of 
ULXs in early-type galaxies and found that "very" lumi- 
nous ULXs (with Lx > 2 x 10 39 erg s _1 ) are consistent 
with the expected number of background sources. This 



result is based in part on Chandra data and uses lumi- 
nosities in the 0.3-10.0 keV bandpass as opposed to our 
choice of the 2-10 keV band. However, as can be seen in 
Figure 12b, even when considering all ULXs in our sam- 
ple (i.e., with Lx > lx 10 39 crg s _1 ), the frequency of 
early-type galaxies containing ULXs is only marginally in 
excess of zero (with ULXs occurring in globular clusters 
as in the case of NGC 1399 contributing to the early- type 
galaxy ULX frequency). 

The main caveat with an analysis of the HRI ob- 
servations is the softness of its bandpass (i.e., 0.1-2.2 
keV). Objects absorbed by column densities in excess of 
~ 10 21 cm -2 are therefore difficult with the HRI. Since 
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columns of this order are typically observed through the 
disks of galaxies, it is likely that ULXs are being missed 
that are close to the nucleus of the galaxies and/or are 
observed in edge-on disks. Chandra and XMM-Newton 
observations will therefore be much more adept at de- 
tecting ULXs in these types of environments. This is 
particularly true for Chandra whose much smaller PSF 
relative to the HRI and XMM-Newton allows nuclear 
ULX to be discerned fro m AGN at sm aller offsets (i.e., 
1" compared with 10"). iSipiorl ()2003[) finds that 8/41 
~ 20% of nearby galaxies in a relatively unbiased Chan- 
dra survey harbor a ULX. This suggests that our HRI 
survey is missing of order ~ 40% of ULXs (albeit with 
large uncertainty due to the size of the Chandra sample) . 
iHornschemeier et al.l 1)2003 ') reports that the frequency of 
ULXs with L x > = 2 x 10 39 erg s'Ms ~ 8% l ocally and 
~ 36% at z ~ 0.1 (|Hornschemeier et al.ll2003j) (suggest- 
ing that the frequency of ULXs may evolve). We are 
similarly finding ~ 8% (see Table ^) which implies that 
any biases in our survey are not particularly large for 
L x >= 2 x 10 39 erg s _1 ULXs, within the errors. 

The lack of any obvious connection between FIR lu- 
minosity (and hence SFR) or the FIR/K band luminos- 
ity ratio distribution and the presence of a ULX sug- 
gests that if such a connection exists, then the biases 
discussed above are particularly acute in star-forming 
galaxies. This may not be surprising given that high-SFR 
galaxies tend to contain significant amounts of molecu- 
lar clouds and dust which would obscure soft X-rays. 
In addition the relatively small size of our sample may 
also be impacting this result. For example, only 23 
galaxies have log Lpm >43.5 erg s _1 and enough HRI 



sensitivity to detect a ULX, and 4 harbor a ULX (af- 
ter "background subtraction" and using a galaxy ra- 
dius of O.5D25). This amounts to a larger ULX fre- 
quency of I7~t\ 2 %, although still statistically consistent 
with the full sample result. Recent Chandra surveys of 
the point source luminosity functions in galaxies sug- 
gest that, particularly at high star- formation rates, the 
point source luminosity function normalization is cor- 
related with the SFR (Grimm et al. 2003, Colbert et 
al. 2004), and that the lumi nosity function is flatter 
in spirals iColbert et al J 12 0041) or (similarly) in higher 
SFR galaxies <|Kilgard et al.ll2002j) . Having both a flat- 
ter luminosity function slope and larger normalization 
at higher SFR strongly implies that higher SFR galax- 
ies (i.e., with SFR greater than several to ten M Q yr _1 
based on Colbert et al. 2004) would tend to harbor multi- 
ple ULXs. This is also suggested anecdotally by the large 
numbers o f ULXs found in merging galaxies such as the 
Antennae l)Fabbiano. Zezas. &: Murravl 12001). Accord- 
ingly the distribution of LpjR and Lpm/L^ for galaxies 
with at least one ULX may not be significantly impacted 
by an increase in the total number of ULXs in galaxies 
with large SFR. However the resolution of these issues 
will ultimately require more detailed statistical analysis 
based on Chandra and XMM-Newton wide-area surveys. 
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